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The direct detection of gravitational waves will provide valuable astrophysical information about 

!>• ■ 

many celestial objects. The SCHENBERG has already undergone its first test run. It is expected 



to have its first scientific run soon. In this work a new data analysis approach is presented, called 



On 
O 

method of independent bars, which can be used with SCHENBERG's data . We test this method 



through the simulation of the detection of gravitational waves. With this method we find the 
source's direction without the need to have all six transducers operational. Also we show that the 
method is a generalization of another one, already described in the literature, known as the mode 
channels method. 
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I. INTRODUCTION 



The detection of gravitational waves will have many implications in Physics and Astro- 
physics. Besides the confirmation of the general relativity theory, it will allow the investiga- 
tion of several astrophysical phenomena, such as the existence of black holes and the mass 
and abundance of neutron stars, thus allowing new scientific frontiers. 

Joseph Weber was the first one to propose feasible gravitational wave detectors in the 
1960's[l]. The antenna he proposed had a cylindrical shape, operated at room temperature 
and was acoustically isolated. Since his pioneering work the detectors based on resonant- 
mass antennas have improved significantly. For the latest generation of such detectors, the 
antenna has a spherical shape. 

The transducer distribution on spherical detectors follows the truncated icosahedron con- 
figuration, proposed by Johnson and Merkowitz in 1993|2] and justified formally by Mag- 
alhaes et al. in 1997[3|. Spherical antennas have omnidirectional sensitivity and all observ- 
ables needed for gravitational astronomy can be obtained from only one spherical detector 
appropriately equipped with six transducers jj]. When fully operational, such detector will 
be able not only to acknowledge the presence of a gravitational wave within its bandwidth: 
it will be able to inform the direction of a source in the sky as well as the amplitudes of 
the wave's polarization components. This is the case of SCHENBERG, the second spherical 
detector ever built in the world and the first equipped with a set of parametric transducers, 
which is installed at the Physics Institute of the University of Sao Paulo (at Sao Paulo city, 
Brazil). It has undergone its first test run in September 8, 2006, with three transducers op- 
erational. Recent information on the present status of this detector can be found in reference 

0. 

There are mathematical models for this detector for the case that all six transducers 
are operational. Such models investigated two situations: one in which the transducers are 
perfectly uncoupled 0, [f], Q] and another in which the transducers are somehow coupled to 
each other 8J. In both cases it is possible to solve the inverse problem and retrieve the four 
astrophysical parameters of interest: the angles that determine the source's position in the 
sky, denoted by (3 (zenith angle) and 7 (azimuthal angle), and the amplitudes of the two 
polarization states as functions of time, h + (t) and h x (t). 

In this work we present a methodology for the solution of the inverse problem in spherical 
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detectors, which we named "the method of independent bars" . Unlike the previous methods, 
it allows the determination of astrophysical parameters using less than six transducers. In 
fact, we will show that it generalizes a previous approach, known as the mode channels 
method. We tested our method with simulations and found encouraging results. In the next 
sections we will show these results in detail. 



II. THE METHOD OF INDEPENDENT BARS 



The determination of directions is sky demands the use of appropriate reference frames. 
For this reason we shall define two frames: the first one is the laboratory system with axes 
defined by the unit vectors (e x ,e y ,e z ). The second frame is the gravitational wave system, 
which will have its axes defined by (ex,ey,e^) with the wave travelling in the ez direction 
and the axes of the polarization amplitudes being (ex,ey)- 

The projection of the polarization state amplitudes of the gravitational signal on the shaft 
of a detector can be expressed by the contraction of two symmetric, trace-free tensors^: 

R(t) = \w(t)D ir (1) 

The response amplitude of one transducer coupled to the spherical detector is given by R(t), 
while h % i(t) is the spatial part of the gravitational wave tensor and is the transducer's 
tensor, which informs the orientation of the transducer on the spherical detector. 

In Appendix |B] we show that the contraction in equation ([I]) is a representation of the 
time-dependent scalar potential given byjf]]: 

1 

<f>{x,t) = -pxjh 3k {t)x k , (2) 
where p is the antenna's mass density and X{ IS cl coordinate location. 



A. The Transducer's Tensor 



It is possible to model the system as if we had one independent bar behind each transducer 



on the detector 



4| . Such a bar detector, whose axis is in the n direction, can be characterized 



by the symmetric tensor 



A 



n w n 



(3) 
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The transducers are located on the antenna and the position of the L transducer is 
described by the unit vector n. In the laboratory system (lab frame) this vector is given by 



riL = sin 9l cos 4>l&x + sin#£, sin 4>L^ y + cos 9ie z , 



(4) 



where 6l and (pL are angles that locate the transducer on the spherical detector. Using 
equations ([3]) and (J4]) the L th transducer's tensor assume the following form in the lab 
frame: 

sin 2 9l cos 2 <Pl \ sin 2 Ql sin 2<j)i \ sin20£, cos</>l 
D L = | sin 2 9 L sin 2(j) L sin 2 6 L sin 2 <p L | sin 29l sin (f) L ■ (5) 
! sin cos 0l | sin 2Q L sin </> L cos 2 9 L 



B. The Wave's Tensor 
Consider the complex null vector m defined by [lOf ] : 

m = — (e x + ze y ) 



(6) 



where and ey are unit vectors in the wave frame. The gravitational wave tensor h can 
be written in terms of these vectors: 



h = 2/i + K(m <g> m) + 2/i x S(m <g) m) 



(7) 



or 



where 



h = h + e + + h x e x , 



(8) 



ex ® e x - e y <g> ey, 
ex <£> ey + ey <g> e x . 



Starting from the lab frame (e^, e y , e z ) it is possible to rotate by 7 around e z to obtain 
(e x /, ej,/, e z r) and rotate by /3 around e y i (Euler ^/-convention) to obtain (e^', e y n, e z n) so that 
e z ii points in the direction of the propagation of the wave e^. 

In the lab frame the vector m thus takes the following form: 

1 



m = —= [(cos (3 cos 7 — i sin j)e x + (cos (3 sin 7 + % cos 7)e y — sin (3e z 
V2 



(9) 



where (3, 7 are the zenithal and azimuthal angles, respectively, that define the direction of 
arrival of the gravitational signal. 

We have determined the expressions for the wave's and the detector's tensors relative to 
the laboratory system. We can now solve the inverse problem. 



C. The inverse problem 



Equation (CQ), valid for one transducer, can be expanded as [HI] 



- [D^h u + 2D[ 2 h 12 + 2D^h lz + D L 22 h 22 + 2^23 + 2D^h m ] = R L 



(10) 



In the case that we have six operating transducers, using the six independent components 
of the spatial part of the gravitational waves tensor, the hij can be obtained as the solutions 
of the system of equations 

'D\, 2D\ 2 2D\, D\ 2 2D\ 3 D\~ 
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(11) 



_D n 2D\ 2 2D\ 3 D\ 2 2D\ 3 D% 3 _ 

This system allow us to calculate the components of the wave's tensor and to check the 
transversality and tracelessness of the wave tensor, which will be present if GR is the correct 
theory of gravitation. On the other hand, in the case the detector is coupled to a set of 
only five transducers (L = 1, ... 5), the last line will be absent. So in its place we use the 
trace- free condition [llj: 

h n + h 22 + h 33 = 0. (12) 

When only four transducers are operational it is still possible to solve the system of equa- 
tions by adding to it the transversality condition of the gravitational wave. This condition 
implies that the spatial part of the gravitational wave tensor is orthogonal to the direction of 
propagation, i.e., h^n^ = 0. In other words, the matrix has one eigenvalue equal to zero 
and, as a consequence, null determinant. We then use the condition of null determinant, 
given by 

hnh 22 h 33 + 2h 12 h 13 h 23 - h 33 h\ 2 - h 22 h\ 3 - h n hj 3 = 0, (13) 



instead of the line that is missing in the system of equations (TTTT) . 

Since we take General Relativity for granted the h matrix, without noise, must be "TT" : 
transverse, equation (TIB"]) , and trace free, equation (fi"2|) . In this case any set of 6 h^ u , 
obtained as described above, should obey such condition. The presence of noise is expected 
to make the h^ u vary so that the transversality and tracelessness of h is not perfect. Of 
course, if the noise is too high it will be difficult to identify the "TT" properties in the h 
matrix. Our ability to extract the signal from the noise will depend on the signal-to-noise 
ratio achieved for the particular event. This has been investigated by Merkowitz 12| and 
Merkowitz, Lobo and Serrano [13], partly using an approach first presented in Magalhaes et 
al. [ljj]. In the present work we will not approach this issue in depth, which will be left to 
be discussed elsewhere. 



III. THE ASTROPHYSICAL PARAMETERS 



With the h obtained from the the system of equations (TTlT) is possible to calculate all the 
relevant parameters of the gravitational wave: the angles (3 and 7 that indicate the direction 
of the wave and the amplitudes of the two polarization states of the wave as functions of 
time h + {t) and h x (t). 



A. The direction of the gravitational wave 

For a gravitational wave that propagates in the direction of the ez axis of its proper 
system of coordinates, in the case of general relativity the transversality condition relation 
holds: 0,Q 

h.e z = 0. (14) 

The wave frame can be written in terms of the laboratory frame through a rotation. 
Following the y-convention for the Euler angles, ez assumes the form: 

ez = sin (3 cos 7 e x + sin (3 sin 7 e y + cos (3 e z . (15) 

Substituting equation (fl5l) in ffT4|) and using the symmetry of we obtain the following 
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system of equations 

/in sin j3 cos 7 + hi 2 sin f3 sin 7 + /ii 3 cos 7 = 0, (16a) 

/112 sin j3 cos 7 + h 22 sin /3 sin 7 + h 23 cos 7 = 0, (16b) 

hi 3 sin cos 7 + h 23 sin /3 sin 7 + /i 33 cos 7 = 0. (16c) 

Since det h = 0, in this system the equations are linearly dependent. Two of them are used 
to obtain (3 and 7: 

1 = Y~h h~h~> (17) 

tan/3=± Yl 3 ~^f 22 V / l + tan^7- (18) 

These two solutions for (3 are related to the two diametrically opposed positions in the sky 
from which the gravitational wave can arrive. 

With the above equations we compute the direction of the incoming wave. It results in 
the same value for each time sample, in the case of a pure (noiseless) signal. In the case of 
a noisy signal, a distribution around the mean values of the j3 and 7 angles is found, as we 
can see from Figures [2] and [3j 



B. Amplitudes of the polarization states 

Using the fact that e + . e + = 2, e x . e x =2 and e + . e x = we obtain h + and h x with 
the aid of equation (jSJ), yielding 

h + =(cos 2 j3 cos 2 7 — sin 2 j)h u + sin 27(1 + cos 2 /3)hi 2 — sin 2/3 cos jh 13 + 

(cos 2 (3 sin 2 7 — cos 2 7)^22 — sin 2(3 sin 7/123 + sin 2 (3h^ 1 
h x = — cos (3 sin 2jhu + 2 cos (3 cos 2^h\ 2 + 2 sin (3 sin 7/113 

+ cos (3 sin 27/122 — 2 sin (3 cos 7/123- 

Once these amplitudes are known other quantities can be calculated. One of them is the 
polarization angle, through the relation 

tan(2^)) = M> 
h+{t) 
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However, this parameter is actually not relevant for gravitational wave detection purposes 
since it is just an arbitrary angle between the wave x axis and an axis where the wave has 
only the instantaneous plus polarization, h + . 

On the other hand, important information can be obtained from the difference between 
the polarization phases. This can be determined if we analyze the time evolution of the 
instantaneous polarizations registered in each sample (the polarizations are time dependent). 
Our ability to determine this phase difference depends on the signal to noise ratio for each 
particular event sampled. Such time varying analysis is beyond the scope of the present 
work. 



IV. TESTING THE METHOD OF INDEPENDENT BARS 



With the objective of testing the efficiency of the method we produced simulated data 
of the SCHENBERG detector. For the waveform we chose a template from the Laboratory 
for High Energy Astrophysics (LHEA) of the NASA/GFSC 16j, one of the three waveforms 



used by Duez and collaborators 



171 . Il8j | and depicted in Figure [p 



This signal was introduced in a simulator of SCHENBERG [191] . that returns the response 
of the six two-mode transducers positioned on the detector their positions obeying the 
truncated icosahedral configuration 20] . The simulation of the detector is possible because 
its mathematical model is known as well as the mathematical models of its noise sources 
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Table [J shows the angles relative to the position of each transducer on the spherical 
detector. We simulated the transducers' responses when detector was excited by a wave of 
the kind shown in Figure [1] arriving at the detector in the direction defined by f3 = 30° and 
7 = 75°. 

In the SCHENBERG simulator the amplitude signal-to-noise ratio (SNR) was estimated 
from Thorne as [22| 

SNU^jy^^lff, (!9) 

where W(f) represents the applied filter function (chosen, in the present work, as the func- 
tion of a pass-band filter). S^if) is the equivalent noise profile (Gaussian, in this case); it 
can be obtained from the detector's output when no useful signal is present, representing 



S 



TABLE I: The angle Ql is the one that a line L, that passes through the center of a pentagonal 
face and the origin of the lab frame, makes with the local zenith and 4>l is the angle that L makes 
with the x axis of the lab frame. 
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79, 1877° 79, 1877° 


79, 1877° 


37, 3774° 


37, 3774° 


37, 3774° 




120° 240° 


0° 


300° 


180° 


60° 



the total noise energy in the sphere. 

The signal-to-noise used was SNR ~ 8 corresponding to a source with behavior as in 
Figure [H located approximately 97 kpc away. It is from this source model that h GW (f) can 
be obtained. 

We used the output of 5 or 6 transducers to test our model, as follows. 



A. Case 1: pure signal 

Once we had the simulated response of the transducers we calculated the parameters of 
the gravitational wave using the method of independent bars. 

The result obtained with the method for the source's direction was exactly the same 
initially imposed: (3 = 30° and 7 = 75°. This shows that the method has good accuracy and 
does not present intrinsic errors. 



B. Case 2: signal plus noise 

We then tested the method in the presence of noise adding Gaussian noise to the data. 
To illustrate the results we projected the sphere onto a plane using the Hammer-Aitoff 



projection 23J], where the direction is given by 9 = 90° — (3 and <fi = —7. 

Figures [2] and [3] show the source's direction as found using the method of independent 
bars for the cases where the data of six or five transducers are used. In those figures the two 
regions of the sky are marked because there is a natural ambiguity in the determination of 
the direction by only one detector. If more than one detector is used then there is a time 
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delay between detections and one can determine whether the wave came from up or down. 

As one can verify, the smaller the the number of transducers' data used in the analysis the 
smaller the precision of the results. In those figures we observe that, for the case in which 
the data of only five transducers are used, the area of interest is larger than in the case 
that use the data of all six transducers. This area is equivalent to the error. The numerical 
results obtained for the method are presented in Table [III 

TABLE II: Numerical results obtained with the method of independent bars for the wave's direction 
in the case of signal in the presence of noise (SNR ~ 8) using five or six transducers' data, n is 
the number of the transducers. 



n 


P 


7 


5 


(27.2 ±7.1)° 


(78.7 ± 5.3)° 


6 


(28.1 ±4.2)° 


(75.3 ±2.9)° 



The values presented in this Table were obtained simulating the wave's incidence 50 times 
and then calculating the average and the standard deviation of the results. 

We have verified the efficiency of the method of independent bars in determining a source's 
direction. Another important result is to show the equivalence between the method of 

n 

independent bars and mode channels model [6J. This will be done in the next section. 

V. EQUIVALENCE BETWEEN THE METHOD OF INDEPENDENT BARS AND 
THE MODE CHANNELS MODEL 

In order to show the equivalence between the mathod of independent bars and mode 
channels model, we must write the transducer's tensor (eqs. OandE]) in its trace-free form: 

Aj = (n® n)ij - -8 i:j , 

where the term |5y ensures the property of tracelessness of the tensor. 

Each component of the matrix D (eq. [TT1) is a function of the angles 9l and <pL that 
informs the position of each transducer. When we impose the trace-free condition on the 
tensor Dij, all this functions can be written as a combination of real spherical harmonics of 
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second order. Therefore the matricial system can be written as 



it 
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(20) 
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where /i m are the quadrupolar representation of the gravitational waves. In Appendix IAI we 
shown the details of this calculation. 

When the transducer are positioned according to the truncated icosahedral configuration 
the matrix of the real spherical harmonics Y is the transpose of the model matrix Bg x5 
proposed by Merkowitz and Johnson 
as 



241 ] . In compact notation equation (TSUI) can be written 

(21) 



R =2 



1 IQn 



Using a property of the model matrix 



15 
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we obtain: 



or, in explicit form: 
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(24) 



R 1 
R 2 
R 3 
R A 
R 5 
R 6 

This result, obtained with the independent bars method, coincides with the one found 
with the mode channels model. Therefore, if the transducers are in a truncated icosahedron 
configuration the two models are equivalent. 

On the other hand, it is worth stressing that even if the transducers are not in this special 
configuration the inverse problem can be solved using the independent bars model. 
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VI. SUMMARY 



The main objective of this work was to show the details of the analytical construction 
of the method of independent bars and its efficiency in the determination of the relevant 
parameters of gravitational waves detected by spherical detector, namely: the two angles 
that determine the source's position in the sky and the amplitudes of the two polarization 
states as functions of time. 

As can be seen from the results, the method of independent bars is efficient for the solution 
of the inverse problem. First we tested the method in the noiseless case and we verified that 
it does not present intrinsic errors, being thus an exact method. Then we tested the method 
adding Gaussian noise to the data, as can be seen in Figures [2] and [3J 

A major advantage of the method of independent bars is its generality, since it does 
not depend on the specific transducers' configuration. The fact that D is a square matrix 
facilitates the inversion of the system's equation, the only requirement being that the matrix 
is non-singular. 

Therefore, the usefulness of this method is evident. For instance, in the absence of any 
transducer coupled to the spherical detector the mode channels mo del [(J does not work. This 
happens because the truncated icosahedron configuration is lost. However, the independent 
bars method still works under this condition. We are currently investigating ways to solve 
the inverse problem with three transducers. We have already got some preliminary results 
that will appear in a future paper. 

In this work we have actually shown that the method of the independent bars is a gen- 
eralization of the mode channels model. The matrix D can be written in terms of the real 
spherical harmonics. In this case, if the transducers' positions obey the truncated icosahe- 
dron configuration we verified that D is the model matrix B 5x6 multiplied by a constant. 

Finally, we point out that this method can be used for the solution of the inverse problem 
in a network of bar antennas, as long as a common reference frame is chosen and time delays 
are accounted for. 
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APPENDIX A 



Each component of the matrix (fTTI) is a function of the angles 9 and <fi that inform the 
position of each transducer. The D functions can be written as combinations of the real 
spherical harmonics as follows: 
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Here the Y£ are defined as functions of the spherical harmonics of second order Yf n 251 ] : 
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Equation ( flOl) can be written in terms of the real spherical harmonics 
R 1 
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In a more compact form this reads 
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where the h m are the quadrupolar representation of the gravitational wave: 

hi = (^22) , 



h 2 
hs 
h i 
h 5 



2 

hi2, 
hi3, 

— %3- 



Therefore the matricial system of the equation fflOl) takes the form: 



R 1 
R 2 
R 3 
R" 
R 5 
R 6 



16n 

"15" 



yl yl 

1 \ 2 2 

I I 1 2 

Y? Yi 

1 1 2 2 

y5 v5 

1 1 2 2 

Y* Yi 



yl yl 
1 3 1 A 



Y 1 



Y 2 Y 2 Y} 



V3 V3 
1 2 1 A 



V 4 V 4 V 4 

J 3 I A J 5 

v5 v5 v5 

2 3 I A I 5 

K, 6 F, 6 



hi 

hi 
h 3 
h.i 
h 5 



which is the same used as equation (T2U1) . 



APPENDIX B 

We want to show in this appendix that equation (J2J) can be written in terms of the 
transducer's response. First we notice that the potential of the tidal force density caused 
for the passage of a gravitational wave is 

<K X > *) = $^ ^pxjh jk (t)x k , 

where p is the mass density and X{ is the coordinate location. The term Xi represents a 
vector for the position so it can be written of the following form: 

Xi 

where e is the amplitude of the vector and n« is represent the unit vector of the position 
vector. Replacing Xi in the equation we get 



(x,t) = ^2^e 2 pn j h jk (t)i 



2 " r ■■ j- .//> \-) n ki 

j,k 
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which can be written in the form 

1 d 2 

We have then the contraction of two tensors: the wave's tensor and a directionality 
tensor that characterizes the plan in which the states of polarization of the gravitational 
wave act. We can impose the traceless condition in riiUj adding the term — Therefore 
the equation assumes the form: 



where the term (jijUk — §<^'fc) is the transducer's tensor of the equation ([3]). Therefore we 
can write 

1 d 2 

^'^E/^ (h jk (t)D jk ), 



or in a more summarized form: 



2 r dt 2 



where 1Z is the transducer's response of the equation (JTJ. Using this equation we can write: 

0(x,£) = a/ y-^Vy^ h m (t)Y m . 
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FIG. 1: Waveform of a gravitational wave emitted by a binary neutron star- neutron star system 
near its last stable orbits. The time and amplitude axes are scaled according to the object masses 
through the multiplication by M/1.5M so i ar , where M so i ar is the Sun's mass. 
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FIG. 2: Determination of the wave's direction using the method of independent bars the presence 
of Gaussian noise with SNR ~ 8. The transducers used are numbers 1, 2, 3, 4 and 6 in Table |I]). 
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FIG. 3: Determination of the wave's direction uning the method of independent bars in the presence 
of Gaussian noise with SNR ~ 8. The six transducers are positioned according to the truncated 
icosahedro configuration shown in Table HI 
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